Electronic transport in an anisotropic Sierpinski gasket 
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We present exact results on certain electronic properties of an anisotropic Sierpinski gasket fractal. 
We use a tight binding Hamiltonian and work within the formalism of a real space renormalization 
group (RSRG) method. The anisotropy is introduced in the values of the nearest neighbor hopping 
integrals. An extensive numerical examination of the two terminal transmission spectrum and the 
flow of the hopping integrals under the RSRG iterations strongly suggest that an anisotropic gasket is 
more conducting than its isotropic counter part and that, even a minimal anisotropy in the hopping 
integrals generate continuous bands of eigenstates in the spectrum for finite Sierpinski gaskets of 
arbitrarily large size. We also discuss the effect of a magnetic field threading the planar gasket on 
its transport properties and calculate the persistent current in the system. The sensitivity of the 
persistent current on the anisotropy and on the band filling is also discussed. 
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I. INTRODUCTION 

Electronic properties of fractal lattices have been 
extensively investigated in the past [3-[2l|, and a 
wealth of knowledge has been accumulated. Such 
deterministic lattices provide excellent examples of 
systems in between perfect periodic order and com- 
plete randomness and exhibit electronic properties 
drastically different from a crystal and a disordered 
material. The variable environment around any site 
of a fractal structure makes such systems quite sim- 
ilar to the random networks. This has motivated 
a series of experiments to examine the fractal be- 
havior, magneto-resistance and the superconductor- 
normal phase boundaries on Sierpinski gasket wire 
networks [^-0 US [HI • The interface of theory and 
experiments has further been strengthened by the 
recent development of non-dendritic, perfectly self- 
similar fractal macromolecules of hexagonal Sierpin- 
ski gaskets (2^ . 

Fractals, in general, exhibit a Cantor set like en- 
ergy spectrum [l| that is highly degenerate, and the 
density of states displays a wide variety of singulari- 
ties. The spectrum may however, broaden up in the 
presence of a magnetic field [sj . The electronic con- 
ductance typically exhibits scaling with a multifrac- 
tal distribution of the exponents One useful 

tool in investigating regular, finitely ramified frac- 
tal lattices has been the real space renormalization 
group (RSRG) methods. The RSRG approach has 
been successfully used over the years to unravel, for 
example, the spectral features of Koch fractals with 
long range interactions Q, localization aspects and 
the length scaling of corner-to-corner propagation 
in fractal glass and several other interesting fractal 
networks [1, • 



Interestingly, RSRG studies on fractal lattices re- 
veal another remarkable property of such systems, 
viz, the existence of an infinite number of isolated 
extended eigenstates that coexist with an otherwise 
Cantor set spectrum. This is non trivial, as the frac- 
tals as such do not have any long range translational 
order. Such extended states sometimes are traced 
back to the existence of multiple cycle fixed points of 
the Hamiltonian, or have been shown to arise out of a 
local positional correlation, giving rise to a resonant 
tunneling effect in some local atomic clusters p3l - 
[26j . However, no conclusive evidence of the possible 
formation of a continuum of states has been obtained 
until very recently, when Schwalm and Moritz [27t 
have presented exhaustive numerical results to show 
that a family of fractals does possess a continuum 
of such states. Incidentally, such continua were also 
'suspected' in some of the previous studies [23 - [26| . 
where an external magnetic field was shown to play 
a role. 

This, to our mind, keeps the problem alive, and 
we look back into the aspect of observing continu- 
ous distribution of eigenstates in an anisotropic Sier- 
pinski gasket (SPG) fractal. The magnetic flux, as 
discussed in the earlier works [23 - [26| through each 
elementary plaquette of an isotropic SPG breaks 
the time reversal symmetry of the electron hopping 
between nearest neighbors. This in turn, already 
generates a kind of anisotropy in the hopping of 
an electron when it travels along an angular bond, 
compared to when it hops along a horizontal bond 
(Fig. H]). Is anisotropy the factor responsible for 
any continuum then ? This question is the central 
motivation behind our present work. Such an is- 
sue had been addressed earlier by Hood and South- 
ern 28[ using a generating function approach cou- 
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pled with RSRG to come to the conclusion that in 
an anisotropic SPG only certain hierarchical states 
persist. 

A Second motivation is of course the fact that, 
grafting a deterministic fractal geometry on a given 
substrate is quite feasible using the present day ad- 
vanced nanotechnology and lithographic methods. 
Therefore, any spectral peculiarity that might come 
out as a consequence of the anisotropy, and is likely 
to get reflected in the transport behavior, would be 
quite possible to observe. This may throw some light 
into the unaddressed issues related to such determin- 
istic fractal lattices, and will definitely enrich the 
experimental aspects. 

In this communication We examine the trans- 
mission spectrum of finite but arbitrarily large 
anisotropic SPG lattices, with and without any mag- 
netic flux threading the system, using an RSRG dec- 
imation method 24]. In addition to this, we also 
calculate the persistent current in such a network, 
and look for any peculiarities in the behavior of the 
current as well as in the Aharonov-Bohm (AB) oscil- 
lations in the transmission that might come out as 
a result of the anisotropy. To the best of our knowl- 
edge, no results are available for persistent currents 
in an SPG network, specially for the anisotropic 
case, and the impact of the scale invariant, mul- 
tiple loop geometry of an SPG on such properties 
are either very little addressed, or unaddressed at 
all. Therefore, a thorough study of the problem, to- 
gether with the literature that already exists will 
possibly help in obtaining some conclusive result 
about the spectra of regular fractal lattices under 
various conditions. 



A 




B H B 



FIG. 1: A Sierpinski gasket and its renormalization to a 
'dimer' 

We find interesting results. Anisotropy indeed 
signals the formation of bands of eigenstates with 
high transmittivity. In general, it is seen that an 



anisotropic SPG is more conducting compared to its 
isotropic counterpart. We have carefully carried out 
an RSRG analysis of the recursion relations of the 
parameters of the Hamiltonian, and observe certain 
unusual flow pattern. The transmission coefficient 
for arbitrarily large lattices has been worked out. In 
addition to this, we present an in-depth study of the 
persistent current observed in a SPG network, and 
the influence of the degree of anisotropy on the cur- 
rent. The behavior of the persistent current is seen 
to be very sensitive to the degree of anisotropy. 

In what follows, we present the results of our in- 
vestigation. In section II, we present the model and 
the RSRG method. Section II contains the discus- 
sion on the transmission spectrum, while, in section 
IV we present the results of the calculation of the 
persistent current. 



II. THE MODEL AND THE METHOD 

We begin by referring to Fig. [T] The Hamiltonian 
of the network is given by, 

H = ^ e„ct c„ + ie'^ ^ cj,c„ + h.c. (1) 

n <nm> 

where, c„ (c„t) is the annihilation (creation) opera- 
tor at the nth site of the gasket, the nearest neigh- 
bor hopping integral t = tx along the 'horizontal' 
bond and it is equal to ty along any 'angular' bond. 
— 27r(/i/(/)o is the phase acquired by the electron 
in hopping along a side of the elementary triangle, 
as shown in Fig. [TJ The phase factor 9 is positive 
if the electron hops around an elementary triangle 
in the clockwise sense, and negative otherwise. We 
have chosen two kinds of 'on-site' potentials depend- 
ing on the local environment of any given site. They 
are symbolized as e and /i, and the 'edge' sites have 
potentials ca and cb corresponding to the vertices 
A and B. For calculating the end-to-end transmis- 
sion we have fixed two semi-infinite ordered chains as 
leads (Fig.l) at the end sites B. The on-site poten- 
tial at the lead-sites has been taken to be constant 
and equal to eo. The nearest neighbor hopping inte- 
gral in the leads is chosen as ^o- 

To obtain the transmission coefficient, as well as 
to judge the nature of the eigenfunctions we renor- 
malize the SPG using a decimation technique [2^ . 
We first explicitly discuss the anisotropic SPG with- 
out any magnetic field. The recursion relations for 
the on-site potentials and the hopping integrals ( in 
the absence of any flux threading the gasket) are 
given by. 
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where, the matrix elements corresponding to the 
'diatomic molecule' are given by, Mn = {E — 

tB,nf/(totx,n) - tx,n/to, Mi2 = -M21 = -{E - 

^B,n)/tx,n, and M22 = -ta/tx,n- Before ending this 
section we just point out that the same decimation 
technique has been used to obtain the recursion re- 
lations in the presence of a magnetic field. The re- 
cursion relations in the latter case are much more 
complicated compared to the field-free case, and we 
decide not to present them here to save space. The 
role of the magnetic field will however, be discussed 
in appropriate places. The results are now presented 
below. 



Using the above set of equations, we renormalize 
an n-th generation SPG (n — f )- times to reduce it to 
a simple triangle. The top vertex is then decimated 
to generate an effective diatomic molecule B — B 
clamped between the leads, as shown in Fig. 1(b). 
The transmission coefficient is then obtained by the 
standard formula [29j . 



(4) 
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III. RESULTS AND DISCUSSION 

A. The transmission spectrum and the 
electron states 

The transmission coefficients for the zero- field case 
and for flux (f) = c/jq/A threading an elementary trian- 
gle have been shown in Fig. [5] for a sixth generation 
SPG in the isotropic limit for the sake of complete- 
ness. As is already known [IJ], the magnetic field 
makes an isotropic gasket more conducting with the 
signature of continuous zones of transmission. When 
a slight anisotropy is introduced, and the magnetic 
field is taken off, it is interesting to observe that. 
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FIG. 2: Transmission across a sixth generation isotropic 
gasket for (a) = 0, and for (b) (j) ~ </'o/4. We have 
chosen ej = for all the vertices of the SPG as well as for 
the sites in the leads, and tx = ty = to — 1 throughout. 



an SPG. In a recent work Schwalm and Moritz [27| 
have discussed precisely this issue in the case of a 
different class of hierarchical lattices. This latest 
observation gives us confidence in regards of the ex- 
istence of bands in fractal lattices which can coexist 
with even a fragmented Cantor like spectrum. 
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first of all the SPG becomes much more conduct- 
ing than its isotropic counterpart, and second, once 
again patches of continuum appear in the transmis- 
sion spectrum. One such example is presented in 
Fig. [3] with tx = 0.9 and ty — 1. With increas- 
ing degree of anisotropy the continua become much 
more pronounced at places. We have not been able 
to provide an analytical proof as to whether real 
continua appear even in such non-translationally in- 
variant systems, but have carried out extensive nu- 
merical investigations involving very fine scan of the 
'apparently' continuous bands found in the transmis- 
sion spectrum. In every case we find the fine struc- 
ture in these 'special' zones free from any real gaps. 
Even, the apparently sharp peaks in the transmis- 
sion spectrum show up finite width when the energy 
is scanned minutely in their neighborhood. Thus, it 
is tempting to claim that anisotropy introduces con- 
tinuous bands of high transmission in fractals such as 



FIG. 3: (a) Transmission across a sixth generation 
anisotropic gasket for <j!> = and with tx ~ 0.9 and 
ty = 1. (b) Fine scan of the energy interval [—1.92, —1.8] 
with SE = 10^7. Other parameters are the same as in 
Fig. 2, and energy is measured in unit of the hopping in 
the leads. 

Are these conducting states extended ? To un- 
derstand this, we have critically examined the flow 
of the hopping integrals tx and ty under successive 
RSRG steps. We find two different types of behav- 
ior. In general, for an arbitrarily chosen energy both 
tx and ty flow to zero quickly under iteration. This 
is understandable, as the spectrum of such a fractal 
lattice is fragmented in general. So, any arbitrarily 
chosen energy, in all probability, would either fall in 
a gap, or would correspond to a localized eigenstate. 
However, Hood and Southern [2^ have eliminated 
the possibility of exponentially localized states in an 
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anisotropic SPG. The flow of the hopping integrals 
under RSRG changes its pattern when we choose an 
energy from what apparently looks like a continuum 
in the transmission spectrum in Fig. [Sfa). For any 
energy within the continuous portion of the spec- 
trum tx remains non-zero for arbitrary number of 
RSRG loops while the hopping ty ultimately flows to 
zero. This feature is always true no matter whether 
the initial value of is less or greater than the initial 
value of ty , and persists for any degree of anisotropy. 
That is, the anisotropy somehow develops a 'prefer- 
ence' for tx, and breaks the symmetry at all scales 
of length. Thus, in the infinite lattice limit the en- 
tire SPG fractal consists of isolated dimers coupled 
by the hopping integral ix,n->oo. For a finite (but 
arbitrarily large) SPG we have coupled the semi- 
infinite leads to the two base atoms (B — B). Thus, 
as ty becomes (practically) zero under RSRG oper- 
ation, we are left with a perfectly ordered 1-d chain 
of atoms with an effective, energy dependent poten- 
tial connected to its nearest neighbor by a hopping 
amplitude tx^n^oo (Fig. S]). So, any energy E that 
will map the original SPG to such a configuration 
will definitely correspond to a conducting state pro- 
vided, the energy falls within the allowed band of 
the leads. What is peculiar about such energy val- 
ues is that, they apparently form a completely con- 
tinuous band. We have scanned the spectrum very 
carefully. One such result of scanning is exhibited in 
Fig. Etb), where we have chosen, quite arbitrarily, 
a certain portion of the spectrum shown in Fig. 2a, 
and scanned the selected zone in an energy interval 
of 10~^. It is obvious that the spectrum retains its 
continuous character. 

Before we end this sub-section, it is worth men- 
tioning that, getting a non-zero transmission is def- 
initely connected to the fact that we have attached 
the leads at the two base atoms of the SPG. Talking 
about an infinite SPG, we should appreciate that we 
can build an infinite SPG using a top down approach, 
when one takes a triangle and keeps on piercing it in 
the SPG design an infinite number of times. We can 
keep the leads attached to the base pair of B atoms 
from the beginning. So, the effective ordered chain 
of atoms which the SPG boils down to, when we 
pick up an energy from the continuum, will always 
be found clamped between the leads. The corre- 
sponding state will be extended in the sense that a 
finite number of sites of the parent SPG will have 
non-zero amplitudes of the wave function. In the 
reverse approach, when one builds an infinite SPG 
by placing SPG's of previous generations on top of 
each other, following the growth rule, it is of course 
not very meaningful to talk about the base atoms 
such as the B — B pair, or even fixing the leads at 



the ends. So, at most the infinite system can break 
up into diatomic clusters as depicted in Fig. 21 
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FIG. 4: Separation of a finite Sierpinski gasket under 
RSRG into detached diatomic clusters and linear ordered 
chain clamped between the leads. 



B. The Aharonov-Bohm oscillations 

The anisotropic SPG displays a wide variety of 
AB-oscillations which are sensitive to the energy 
chosen, as well as on the degree of anisotropy. The 
period of oscillations has been found to be equal to 
(j>o/2 if we choose, say, E — 0. The AB-oscillations 
are displayed in Fig. [5] for the cases {tx — ty — 1), 
{tx = l,ty = 2) and {tx — l,ty = 10). The simple 
oscillation profile in the isotropic case develops into 
a profile with multiple peaks and valleys (Fig. [5jb)) 
as the anisotropy grows larger and finally, in the case 
of very large degree of anisotropy, the AB-oscillation 
profile consists of delta-like peaks at special values 
of the magnetic flux, indicating that, the anisotropic 
gasket triggers ballistic transmission for E = only 
at certain special values of the magnetic flux. 



C. The persistent current 

It is well known that a simple Aharonov-Bohm 
flux threading a metallic or semiconducting ring gen- 
erates an equilibrium persistent current [s^]- 33 1. 
With reference to our SPG fractal, where the same 
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which shows multiple band crossings, as observed 
by Rammal and Toulose [2] earlier. The band cross- 
ings, related to the symmetry and degeneracy of the 
eigenvalues at various values oi (j)/4>o are already dis- 
cussed in details in reference 
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FIG. 5: (a) Aharonov-Bohm oscillations in a 4th gener- 
ation SPG for E = 0. (a) Isotropic SPG with — I and 
ty — 1, (b) anisotropic SPG with tx = 1, and ty — 2, 
and (c) anisotropic SPG with tx = 1, and ty — 10. The 
potential and other parameters are the same as in the 
previous figures. 



flux penetrates each elementary triangular plaque- 
tte, the eigenvalues and the eigenstates are flux pe- 
riodic with period </)o,and hence, as is well known, 
the equilibrium persistent current carried by the en- 
ergy level En is given by. 



I71 = -c 



(5) 



and, the total current in the ring is, = J2n=i 
where, Ng is the total number of electrons in the 
system. It is clear that the nature of the energy 
spectrum plays a major role in dictating the persis- 
tent current in a system. We therefore first have a 
look at the flux dependence of the eigenvalue spec- 
trum of a 4-th generation SPG fractal (Fig. [S]). For 
the sake of comparison, we present the spectrum of 
the SPG in the isotropic limit (Fig. |6l top panel). 
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FIG. 6: Flux dependence of the eigenvalue spectrum of 
a 4th generation Sierpinski gasket, (a) The isotropic 
gasket with tx = ty = 1, (b) the anisotropic gasket with 
tx — 1, ty = 2, and (c) with tx = 1 and ty — 10. In each 
cast the on-site potential is taken to be zero throughout. 



With the introduction of anisotropy in the hop- 
ping integrals { ^ ty ) gaps start to open up in 
the spectrum lifting the degeneracy in places. As 
the degree of anisotropy is increased, the spectrum 
shows clear signature of having groups of states ly- 
ing close in energy and separated by wider global 
gaps, that bring a flavor of sub-band structures in 
the spectrum. The middle and the bottom panels 
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FIG. 7: Persistent current in a 4th generation Sierpinski 
gasket, (a) The isotropic gasket with tx = ty = 1. The 
solid and the dashed curves correspond to A*'e = 4 and 
15 respectively, (b) The anisotropic gasket with = 1, 
ty — 2 (solid line), and tx = 1, ty = 10 (dashed line) 
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(dashed line). In each cast the on-site potential is taken 
to be zero throughout. 



The flattening of the bands in this case is notewor- 
thy, which finally brings in non trivial changes in the 
persistent current. 

We have examined in details the variation of 
the persistent current for both the isotropic and 
anisotropic SPG's. for the isotropic case with = 
ty = I, the current does not exhibit serious devia- 
tions (in its qualitative features) from that in a sin- 
gle loop ^ as long as is small. With increas- 
ing values of Ne, a slight rounding off of the I{4>) 
is observed. The overall current keeps on increas- 
ing as A^e increases. The first variation shows up for 
A^e = 10, when additional kinks appear at 4> = </'o/2. 
With increasing A^e the SPG starts displaying fea- 
tures different from what one finds in simple loop 
structures. The features persist as we approach the 
half filled case which, in this study is A'e = 21, when 
a precise 0o/2 periodicity in /(0) is observed. In 
Fig-EIa) we show the variation of the persistent cur- 
rent when A^e = 4 (solid line) and A^g = 15 (dashed 
line) . 

Introduction of anisotropy in the values of the 
hopping integral brings in different features in the 
/((/))-0 curves. Though several basic features such 
as the ipo periodicity, or the special case of 0o/2 pe- 
riod at Ne = 21 persist, the degree of anisotropy, 
i.e. the relative values of and ty strongly influ- 
ence the current. Needless to say, the filling fac- 
tor also plays a crucial role in determining the cur- 
rent profiles. For example, we have found that with 
(tx,ty) — (1, 10), the current drops compared to the 
case when {tx,ty) = (1,2) when we take A^e = 21. 
This has been tested with several values of the hop- 
ping and we come to the conclusion that at half fill- 
ing, the effect of anisotropy is to reduce the overall 
current. The current profile becomes smooth with 
the discontinuous features of isotropic cases for low 
A^e as we introduce anisotropy. This is mainly due 
to the gaps which open up in the spectrum as the 
anisotropy is introduced. In this aspect, anisotropy 
plays a role equivalent to the role that disorder plays 
in one dimensional rings. 
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